Pulling gold nanowires with a hydrogen clamp 
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In this paper an experimental study of the interaction of hydrogen molecules with gold nanowires 
is presented. Our results show, that chains of Au atoms can also be pulled in hydrogen environment, 
however in this case the conductance of the chain is strongly reduced compared to the perfect trans- 
mission of pure Au chains. The comparison of our experiments with recent theoretical prediction for 
the hydrogen welding of Au nanowires implies that a hydrogen molecule can even be incorporated 
in the gold nanocontact, and this hydrogen clamp is strong enough to pull a chain of gold atoms. 

PACS numbers: 73.63.Rt, 73.23.-b, 85.65.+h 
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I. INTRODUCTION 

The behavior of gold nanojunctions have been widely 
investigated in the past years. Due to its inertness, 
gold is an easy-to-use material, and Au nanojunctions 
can even be studied under ambient conditions with a 
low risk for contamination. This monovalent metal was 
found to be a perfect candidate for studying the quan- 
tum nature of conductance in a single- atom contact. It 
was shown by various experimental and theoretical meth- 
ods that a monoatomic Au junction has a single conduc- 
tance channel with perfect transmission. 1 Surprising was 
the discovery that Au has a special property of form- 
ing monoatomic chains during the elongation of a single- 
atom contact^ 

Recently, another remarkable property of gold was pro- 
posed: various theoretical simulations have pointed out 
that - contrary to the well-known inertness of macro- 
scopic Au surfaces - gold nano structures are strongly re- 
active constructions of nature This proposal is sup- 
ported by experimental studies on the catalytic proper- 
ties of Au nanoclusters£ The possibility for the chemical 
interaction of gold nanowires with adsorbants was ini- 
tially advised in order to explain the anomalously large 
inter-gold distances in electron microscope images of Au 
nanochains In this paper we present an experimen- 
tal study on the chemical interaction of gold nanochains 
with the simplest molecule, H2. We compare our results 
to recent theoretical calculations investigating the possi- 
ble configurations, where hydrogen is chemically bound 
to the Au nanowire^ 

Our measurements were performed on high-purity gold 
samples with the mechanically controllable break junc- 
tion (MCBJ) technique under cryogenic circumstances. 
We have investigated the characteristic behavior of the 
nanocontacts by acquiring a large number of conductance 
traces during the repeated separation of the electrodes, 
and collecting the data into conductance histograms^ A 
typical conductance histogram measured in ultra high 
vacuum is presented in Fig. Ui. The most pronounced 
feature is a sharp peak at the quantum conductance unit, 
Go = 2e 2 /h. This peak arises due to flat conductance 
plateaus at IGq corresponding to single atom Au contacts 
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FIG. 1: Panel (a) is the conductance histogram of pure gold 
in high vacuum, exhibiting a sharp peak in the IGo region. 
This peak arises due to flat conductance plateaus at IGo corre- 
sponding to single atom Au contacts or monoatomic Au wires 
(the latter is shown in the inset, an interatomic Au-Au dis- 
tance in the chain is 2.5 A). The area graph in Panel (b) is 
the histogram of gold in hydrogen environment, showing the 
appearance of new atomic configurations in the low-G region. 
The inset shows an example for a trace with a well- structured 
periodic behavior in the low- conductance region, which is at- 
tributed to the pulling of a gold chain with an incorporated H2 
molecule. The line graph in panel (b) is the histogram for the 
hydrogen affected traces (see text). 



or monoatomic Au wires (see inset) both having a single 
conductance channel with almost perfect transmission. 1 
The lower edge of the histogram shows that nanocontacts 
with lower conductance can also occur down to 0.7 Go, 
but below that value the histogram has no counts. In 
hydrogen environment, however, a low-conductance tail 
arises in the histogram [area graph in Fig. Us], that is 
the conductance of the nanojunctions can take any value 
below the quantum unit. The growth of the tail in the 
histogram is attributed to the appearance of new, low- 
conductance atomic arrangements in the hydrogen envi- 
ronment. Usually a peak near 0.5 Go is also superposi- 
tioned on the low- conductance tail with varying ampli- 
tude, which was in the focus of our previous workiiS The 
object of this study is to understand the microscopic na- 
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ture of the new hydrogen-related atomic configurations 
through the statistical analysis of the conductance traces. 
The emphasis is put on the whole low-conductance tail, 
and the results are insensitive to the weight of the peak at 
0.5 Go- First we show with the investigation of plateau's 
length histograms that gold nanocontacts in hydrogen 
environment are generally exhibiting the formation of 
atomic gold chains, however, the interaction with hydro- 
gen can strongly reduce the transmission of the chain. 
Next, we demonstrate a unique feature, which is ob- 
served on a portion of the conductance traces showing a 
well-structured periodic behavior. An example for such 
traces is already presented in the inset of Fig.^D. Based 
on recent theoretical simulations we attribute this phe- 
nomenon to the pulling of gold nanowircs with an incor- 
porated hydrogen molecule. 

Our measurements are in close relation to recent stud- 
ies on the interaction of hydrogen molecules with plat- 
inum and palladium nanocontactsii*i£ In both cases the 
histogram in high vacuum shows a pronounced peak at 
~ G = 1.5 Go, which is attributed to the conductance of 
single atom Pt and Pd contacts. The interaction with 
hydrogen, however, completely suppresses this peak, and 
new, hydrogen-related peaks grow. In the case of plat- 
inum the vibrational modes detected by nonlinear I — V 
characteristics gave clear evidence that a molecular hy- 
drogen bridge is formed between the platinum electrodes, 
which has perfect transmission through a single conduc- 
tance channel. In palladium the transmission is reduced 
due to the dissolution of H atoms in the Pd electrodes. 

In gold the novel phenomena due to the interaction 
with hydrogen is reflected by the growth of the low- 
conductance tail in the histogram, but the presence of 
hydrogen does not destroy the sharp peak at the quan- 
tum unit, which is attributed to single atom Au contacts 
or monoatomic Au chains. It implies that in gold the 
interaction with hydrogen can coexist with the charac- 
teristic behavior of pure gold junctions. 



II. CHAIN FORMATION IN HYDROGEN 
ENVIRONMENT 

As it was presented in Fig. ^ hydrogen induces new 
low-conductance configurations. Further on it will be 
demonstrated by the statistical analysis of the conduc- 
tance traces that the observed new configurations are also 
related to atomic gold chains. 

To make a quantitative analysis possible, we define two 
conductance regions, as shown by the shadowed areas in 
Fig. n The first region ("IGo region") is identified with 
the interval where the peak at G = 1 Go grows in the con- 
ductance histogram of pure gold, G/Gq 6 [0.7,1.1]. The 
second, so-called "low-G" region is the interval, where 
the new hydrogen-related configurations arise in the his- 
togram, G/G S [0.05,0.7]. 

We found that even in hydrogen surrounding a lot of 
traces show the behavior of pure gold, exhibiting plateaus 
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FIG. 2: Plateaus' length histograms for gold junctions in hy- 
drogen environment. Panel (a) shows the plateaus' length his- 
togram for the traces that are not affected by hydrogen. Panel 
(b) contains plateaus' length histograms for the hydrogen- 
affected traces (~ 50% of all curves). The line graph shows 
the distribution of the lengths in the IGo region, the dashed 
line shows the distribution of the lengths in the low-G region, 
whereas the area graph presents the distribution of the summed 
lengths of the low-G and IGo regions. 

near 1 Go without any data points in the low-G region. 
The histogram for these curves completely resembles the 
one taken in the absence of hydrogen. Next we investi- 
gate the other portion of the traces that do show low-G 
configurations. We shall call a trace "hydrogen- affected" 
if it stays for a reasonable length of > 0.25 A in the low-G 
region. 

With these definitions we find that typically ~ 50% 
of the traces is hydrogen-affected. The histogram for 
this portion of the traces is plotted by the line graph in 
Fig.OJj. This histo gram can be directly compared to the 
one for the whole data set (area graph), as both of them 
are normalized to the number of traces included. Natu- 
rally, for the hydrogen-affected curves the weight in the 
low-G region is significantly larger. At G > 1.1 Go the 
two histograms precisely fall onto each other, whereas the 
peak near the conductance quantum is significantly sup- 
pressed for the hydrogen-affected traces. It implies, that 
the appearance of the new low-G configurations does not 
affect the evolution of the conductance of larger junc- 
tions, but it strongly influences the behavior of single 
atom contacts or atomic chains. A deeper insight into 
this correlation can be obtained by investigating the dis- 
tribution of the lengths of the traces in the above defined 
two regions. 



A. Plateaus' length analysis 

Figure presents the distribution of the lengths in 
the IGo region for those traces that are not affected by 
the hydrogen. This so-called plateaus' length histogram 
shows an equidistant peak structure, which is the proof 
for the formation of chains of gold atoms, as it was shown 
by Yanson et ali£ The peaks in this histogram correspond 
to the rupture of chains with different number of gold 
atoms included, and the distance between the peaks cor- 
responds to the interatomic Au-Au distance in the chain. 
This plateaus' length histogram is identical to the one 
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FIG. 3: The figure shows the determination of the lengths 
in the different conductance regions for an individual conduc- 
tance trace. The length of a trace is measured by the num- 
ber of data points in a certain region. The summed length 
means the length in the entire IGo and low-G region, i.e. 
G/Go G [0.05, 1.1]. 



taken before admitting hydrogen to the junction, thus 
it can be regarded as a reference histogram showing the 
formation of pure gold atomic chains. The peaks in this 
histogram are also used to calibrate the length scale, as 
in a pure gold atomic chain the interatomic distance was 
measured to be 2.5 

Next we investigate how this chain formation is influ- 
enced by the interaction with hydrogen. In the case of 
pure gold the conductance of the last atomic configura- 
tion before the breakage always shows a flat plateau in 
the close vicinity of IGo (see the inset in Fig.UJi), and ac- 
cordingly the plateaus' length histogram is built from the 
distribution of the lengths in the IGo region. In the case 
of the "hydrogen-affected" curves, however, the conduc- 
tance can take any value below the quantum unit, thus it 
makes sense to study the distribution of the lengths both 
in the IGo and in the low-G region. The determination 
of these lengths is demonstrated in Fig. [21 

The proper analysis of the distribution of the plateau's 
length gives distinctive information about the possible 
microscopic choreographies of the rupture. For instance 
a strong interaction with hydrogen might completely de- 
stroy the chain formation, which would result in the 
disappearance of the equidistant peak structure in the 
plateau's length histogram, and only the first peak would 
survive. As another possibility, one can consider that 
perfectly transmitting chains are formed just like in pure 
environment, and the new low- conductance configura- 
tions are only appearing after the rupture of these chains. 
In this case the distribution of the plateau's length in 
the IGo region should show the same equidistant peak 
structure as the reference histogram in Fig. [2b- In con- 
trast, the measurement shows that for the hydrogen af- 
fected traces the plateau's length histogram for the IGo 
region is strongly distorted compared to the reference 
histogram (solid line in Fig. 12b). A plateau's length his- 
togram can also be plotted for the low-G configurations, 
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FIG. 4: Analysis of the plateaus' length histogram for the 
summed length region. The area graph presents the measured 
distribution of the summed lengths of the low-G and IGo re- 
gions. The dashed line shows the calculated distribution of the 
summed length which means the convolution of the plateaus ' 
length histograms for the low-G region and IGo region assum- 
ing that the lengths in these two regions are independent. The 
appearance of the peaks demonstrates the correlations indicat- 
ing the chain formation in the entire low-G and IGo regions. 



exhibiting a completely featureless distribution (dashed 
line in Fig. |5p). Surprisingly, the distribution of the 
summed length of the low-G and IGo region [area graph 
in Fig. [2b] shows the well-defined peak structure just like 
the reference histogram [Fig. [2b] ■ (For the definition of 
the summed length see Fig. [21) Observing peaks at the 
same equidistant positions as in the case of pure gold 
demonstrates that the chain formation is preserved in 
hydrogen environment. But contrary to pure gold, the 
peaks appear in the plateaus' length histogram of the 
entire IGo an d low-G region which also has an impor- 
tant consequence that the new low-conductance configu- 
rations are not arising after the rupture of the chain, but 
both the IGo and the low-G configurations are a part of 
the chain formation. This can be explained as follows: an 
atomic gold chain is formed just like in pure environment, 
but during this process at a certain point the hydrogen 
binds to the chain, which results in the strong reduction 
of the transmission. An individual trace corresponding 
to this process is shown in Fig. [3 It is noted that the 
interaction with H2 strongly reduces the conductance of 
the chain, but the atomic periodicity agrees with that of 
pure Au chains. 

As the length of the entire IGo and low-G regions fol- 
lows a well-defined, periodic distribution, the lengths in 
its two subregions must be correlated. This correlation 
is demonstrated in Fig. 0] If the length of the configu- 
rations in the low-G and IGo region were independent, 
the distribution of the summed length could be expressed 
simply as the convolution of the distribution functions of 
the two subregions, which is shown by the dashed line. 
The calculated convolution curve assuming the indepen- 
dence of the subregions only gives the rough tendency of 
the measured one (area graph), but it can not explain 
the presence of the peaks. It nicely demonstrates again 
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FIG. 5: The average returning length as a function of the 
pulling length for the hydrogen-affected traces (circles) and 
the non-affected traces (diamonds). The pulling length is 
the length of the plateau in the [0.05, 1.1] region before break, 
whereas the returning length is the length measured from point 
of breakage back to point, where a contact with G > 0.05 Go 
is reestablished. 



that the low-G region and IGo region must be considered 
together and the entire G/Go € [0.05,1.1] is the physi- 
cally relevant conductance region for the chain formation 
process in hydrogen environment. 

Further confirmation for the chain formation process 
for the hydrogen affected curves can be obtained by the 
statistical analysis of the returning lengths following the 
idea of Yanson et al. 2 If an atomic chain is being broken 
the atoms in the chain collapse back to the electrodes, 
thus in order to reestablish the contact the electrodes 
must be pushed towards each other by almost the same 
distance as the length of the chain was. This means that 
the relation between the length of the plateau before the 
break and the so-called returning length is approximately 
1:1. Fig.JSJshows that this relation holds not only for the 
traces that are non-affected by the hydrogen but also for 
the hydrogen affected curves. 



B. New class of chains 

In previous studies solely the analysis of the distribu- 
tion of the plateaus' length was applied to demonstrate 
the chain formation process and to determine the atomic 
periodicity of the chains, as generally the chain forma- 
tion is not visible on the shape of individual conduc- 
tance traces. In hydrogen environment, however, we have 
found that a part of the traces (~ 5% of all curves) ex- 
hibit a well structured periodic behavior in the low-G 
region, as demonstrated by the curves in Fig. |SJ These 
traces show periods with positive slope, where the con- 
ductance is significantly increasing as the electrodes are 
pulled apart. At the end of these periods the conductance 
frequently jumps to the initial value, and the same behav- 
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FIG. 6: Conductance traces showing a well-defined periodic 
structure in the low-G region. The bottom curve is on true 
scale, while the other curves are artificially shifted on the con- 
ductance axis for the sake of visibility (the corresponding zeros 
are signed on the right axis). The interatomic distance in the 
pure gold chain is indicated by grid lines. Approximately 5% 
of the traces shows this special behavior. 
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FIG. 7: The average of all the traces that have a length > 8.8^4 
in the low-G region. The averaging was done as a function of 
the length measured backward from the point of breakage. 



ior is repeated several times. Generally the conductance 
of these periodic sections is situated in the bottom of the 
low-G region, that is in the interval G/Gq 6 [0.05,0.3] 
(which will be called as "tail region" ) . The peak-to-peak 
amplitude of these oscillations can be as large as 0.2Go, 
which frequently exceeds the mean value of the conduc- 
tance. The grid lines of the figure are drawn with a sep- 
aration corresponding to the interatomic distance of the 
pure gold chains, 2.5A. It is apparent that the length of 
these periods is close to the interatomic gold-gold dis- 
tance in a monoatomic chain, which implies that these 
special conductance traces are also related to the chain 
formation. 

To check the statistical relevance of the behavior 
demonstrated by the curves in Fig. [HI we have calculated 
an average conductance curve [Fig.|7j. In order to deter- 
mine the average behavior along several interatomic dis- 
tances, we have restricted the averaging for those traces 
that stay for a length of > 8.8 A in the low-G region. The 
average conductance curve also shows the oscillatory be- 
havior with a period coinciding with the interatomic gold- 
gold distance, showing that the well structured periodic 
behavior of the traces in the low-G region is a substantial 
feature of gold nano junctions breaking in hydrogen envi- 
ronment. We attribute this behavior to a unique type of 
contact breakage, corresponding to a special microscopic 
process, which will be discussed in Sec. IIIII 
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C. Comparison of the two types of chains 

The plateaus' length analysis discussed in Sec. Ill Al 
has shown that gold nanocontacts in hydrogen envi- 
ronment are generally showing the formation of atomic 
gold chains, however the interaction with hydrogen can 
strongly reduce the transmission of the chain. In sec- 
tions III BI Figs. HO and [7| demonstrated a well-structured 
periodic behavior, which is observed only on a portion of 
the traces, and which is also related to some kind of chain 
formation process. In the latter case the chain formation 
is not just statistically detectable, but it is seen on the 
shape of the individual conductance curves. In the fol- 
lowing part it will be analyzed how these two phenomena 
are related to each other. 

It is to be emphasized that not only the periodicity 
but also the pronounced positive slope is an interesting 
peculiar property of the conductance curves presented in 
Fig. El The positive slope of a conductance trace means 
that the transmission of a certain atomic configuration 
is increasing despite of the fact that the electrodes are 
pulled apart and so the average distance between the 
atoms is increasing. Such an unusual behavior is not 
observed in pure gold junctions. Surveying the individ- 
ual conductance traces in hydrogen surrounding one by 
one we found that practically only the curves presenting 
the well-defined periodic structure exhibit positive slope. 
Therefore, this is a distinctive feature, which can be used 
to statistically separate the special periodic traces from 
the rest of the hydrogen affected curves. We found that 
the periodic traces can be selected with good certainty by 
sorting out those curves, for which the length of the parts 
with positive slope in the low-G region is larger than 30% 
of the the whole length in the low-G region. With this 
analysis one can plot the plateaus' length histograms sep- 
arately for the well structured periodic traces, and for the 
rest of the curves, as shown in Fig. El The dashed curves 
are the plateaus length histograms for the well-structured 
periodic traces with pronounced positive slope, whereas 
the solid lines show plateaus length histograms for the 
rest of the hydrogen affected traces. 

Panel (a) in Fig. |H] presents the plateaus' length his- 
tograms for the summed length in the low-G and IGo 
region. The traces without pronounced positive slope 
(solid line) show the same behavior, which was already 
demonstrated in Fig. |2 an atomic chain formation pro- 
cess taking place in the entire conductance region of 
G/Go — [0.05,1.1]. Contrary, the traces with well- 
structured periodic behavior do not show this feature, 
the plateaus length histogram for these curves is com- 
pletely featureless. 

As the special periodic behavior basically appears in 
the tail region, it is interesting to plot the plateaus' length 
histograms for the tail region as well (Fig. |SJd) . In this 
case the distribution of the lengths for the periodic traces 
shows equidistant peaks separated by a distance of 2.5 A 
whereas the distribution of the lengths for the rest of the 
curves is completely featureless. 
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FIG. 8: Plateaus' length analysis for different types of hy- 
drogen affected traces. The plateaus ' length histograms shown 
by the dashed/solid lines are built from the traces that have 
longer/shorter than 30% section with positive plateaus' slope. 
In panel (a) the distribution of the summed length in the en- 
tire IGo and low-G region is plotted, while panel (b) shows the 
distribution of the lengths in the tail region, G/Go £ [0.05, 0.3] 



This analysis shows that two different types of hy- 
drogen related chain formation processes can be distin- 
guished. The first one is seen for the majority of the 
hydrogen affected traces, during the pulling of a chain at 
a random point the conductance switches from the 1 Go 
value to the low-G region, but the entire length remains 
the same as in pure gold. This is demonstrated by the 
solid line in Fig.|SJi. The second one is only observed for 
a portion of the traces (~ 5%), the conductance of the 
chain drops to the tail region where a well-structured pe- 
riodic behavior is visible on every subsequent gold-gold 
interatomic distance (see the dashed curve in Fig. [SJ . In 
the next part the possible microscopic background for 
these two processes is discussed. 



III. TOWARDS THE MICROSCOPIC PICTURE: 
COMPARISON WITH THEORY 

Further steps towards the microscopic understanding 
of the observed phenomena can be done by comparing 
our experimental results with recent computer simula- 
tions on the interaction of hydrogen molecules with gold 
nanocontacts. Barnett et al. have found that the ba- 
sically inert gold surface becomes highly reactive with 
hydrogen at nanoscale due to the relativistic corrections 
in the electronic structure. 4 Various atomic configura- 
tions were studied, including a hydrogen molecule weakly 
bound to the side of a gold chain, a H2 molecule incorpo- 
rated in the gold chain, and dissociated H atoms sitting 
in the chain. This study has shown that a H2 molecule 
incorporated in the gold chain is a robust configuration 
with high binding energy, which can be reached without 
any barrier. All these configurations have conductances 
in the range of — 0.25 Go, which coincides with the inter- 
val, where the strongest weight is observed in the low-G 
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FIG. 9: Illustration for the weakly bound hydrogen molecule 
(a) and the mechanism of pulling an atomic gold chain with 
a hydrogen molecule (b). 



region of the hydrogen-affected conductance histogram 
[Fig. Eft. 

The configuration where the hydrogen molecule is 
weakly bound to the side of the chain (see Fig. 
might give an important contribution to the general be- 
havior demonstrated by the plateau's length histogram 
[Sec. Ill Al Fig. The formation of atomic gold chains is 
preserved, but the chain gets distorted due to the binding 
of the H2 molecule and thus its transmission is reduced. 

On the other hand the configuration, where a hydro- 
gen molecule is incorporated in the chain provides a rea- 
sonable explanation for the pronounced positive slope of 
a part of the conductance traces in our measurements 
[Sec. IHBl Fig. EJ. The simulations have shown that the 
molecule can be oriented in different angles with respect 
to the contact axis. The perpendicular configuration 
has a smaller conductance of ~ 0.1 Go, but as the elec- 
trodes are pulled apart and the molecule turns towards 
the parallel configuration the conductance grows up to 
~ 0.25 Go- During this process the force acting between 
the molecule and the gold nano junction reaches values as 
high as ~ 0.5nN, and before the complete rupture it can 
even grow to ~ 1.4nN^ Note that these forces are in the 
same range as the ones required for pulling monoatomic 
gold chains (0.7 - 1.5nN)i4 

Based on these considerations we give an explicit pro- 
posal for the well-defined periodic behavior of the con- 



ductance traces demonstrated in Fig. H3 We suggest 
that as a gold nano junction is pulled apart, a hydrogen 
molecule gets incorporated in the chain. During further 
retraction of the electrodes this hydrogen clamp does not 
break, but it is strong enough to pull further gold atoms 
into the chain. Once a new gold atom jumps into the 
chain the molecule gets in a tilted position with smaller 
conductance, but as the contact is pulled the molecule 
straightens out and thus the conductance increases until 
the point where the next atom jumps into the chain. An 
illustration of this process is shown in Fig. 

We have also performed I — V curve measurements 
in order to determine the vibrational energies of the 
molecule in the chain similarly to Ref.0. We have found, 
that in the low-G range the nonlincarity due to quantum 
interference phenomenai£*i£ is so high that the vibra- 
tional spectrum could not be resolved. 

IV. CONCLUSIONS 



In conclusion, we have shown that hydrogen molecules 
can strongly interact with gold nano junctions, which 
is reflected by the appearance of new, low-conductance 
atomic configurations. By the analysis of plateaus' length 
histograms we have demonstrated that these low-G con- 
figurations are not established after the formation of 
monoatomic gold chains, but they are a natural part 
of the chain formation process. In addition we have 
shown, that the conductance traces frequently show a 
well-defined periodic behavior with positive slope in the 
low-G region. This observation is suggestive for a spe- 
cific chain formation process: a hydrogen molecule gets 
incorporated in the gold nanojunction, and this hydrogen 
clamp is strong enough to pull a monoatomic gold chain 
from the electrodes. 

This work has been supported by the Hungarian re- 
search funds OTKA F049330, TS04988L A. Halbritter 
is a grantee of the Bolyai Janos Scholarship. 
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